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The nonstandard picture of a turbulent field is presented in the article. By the concepts of non- 
standard mathematics, the picture describes the hierarchical structure of turbulence and shows the 
mechanism of fluctuation appearing in a turbulent field. The uncertainty of measurement is pointed 
out. And the fundamental equations of turbulence are given too. It is natural in this picture that 
the reasonable closure methods can be obtained and seem to be precise. 

1 Some Concepts of Nonstandard Analysis 

The nonstandard analysis theory of turbulence will be presented in this 
paper. The nonstandard analysis means the mathematic fundament of the 
nonstandard picture of turbulence. Therefore, some concepts of nonstandard 
analysis should be introduced in the first place. 

On the basis of mathematic logic, A.Robinson |lj proved in the sixties of 
the last century that a real number system R can be expanded into the hyper- 
real number system R* . Besides all real numbers (i.e., standard numbers), R* 
contains also hyperreal numbers (i.e., nonstandard numbers). The infinitesi- 
mal e and infinite L are elementary nonstandard numbers. There are other 
nonstandard numbers, for example, ± e, ^e, ^L, ± L, ee. LL. etc. stands 
for any real number). 

It is known that real numbers one-to-one correspond to the points on 
a real line, but on the real line there are no such points which correspond 
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to nonstandard numbers. In other words nonstandard numbers can not be 
indicated by the points on a real hne. The concept of monad needs to be 
presented to indicate the nonstandard numbers. For every real number ^, 
there exists a monad that is composed of the numbers infinitely close to ^ and 
^ belongs also to the monad. In fact, there exist infinite monads containing the 
real number ^. The size of all these monads is infinitesimal £. The numbers in 
this monad, except ^, are all nonstandard numbers of which ^ is the standard 
part. 

So any point on a real line is expanded into a monad. The real number, 
represented by the point on the real line, is the standard part of the nonstan- 
dard numbers contained in the monad. And a real line, formed by the points 
which indicate only real numbers, is expanded into a hyperreal line composed 
of monads. The set of the standard parts of all these monads composes the real 
line mentioned above. Every monad on a hyperreal line is also called a point 
(i.e., standard point). Clearly every such standard point is not the absolute 
point which has the length of absolute zero, but the monad whose length is 
infinitesimal e. 

Every monad must contain only one real number. If a monad contains the 
real number ^, we call the monad as ^-monad. 

The physical world is in a hierarchical structure. This hierarchical struc- 
ture can be described by the concepts of nonstandard mathematics. For ex- 
ample, there are three levels: level a (with infinitesimal e, infinite s"^), level /3 
( with infinitesimal infinite e) and level 7 ( with infinitesimal £~^, infinite 
e~^). Any monad on a hyperreal line in level a is also a whole hyperreal line 
in level /3. However, the whole hyperreal line in level a is only a monad of the 
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hyperreal axis in higher level 7. Prom the viewpoint of hierarchical structure, 
the meanings of infinite and infinitesimal are not fixed. It is understandable 
in physics. Physicists, in fact, can even determine how large infinite and in- 
finitesimal arc in a specific physical problem, and often take a certain number 
as infinite or infinitesimal. In physics any number (not absolute zero) is taken 
as infinite, infinitesimal or finite in various specific condition. 

The coordinate systems in two levels will be set up for the convenience 
of description. Let space coordinates be {xi,X2,X3), time t in level a and 
space {x[,X2,x'^), time t' in level /3, the range of values of {x[,X2,x'^,t') be 
[0, Li), [0, L2), [0, L3) and [0,T) respectively. Here Li,L2,I/3 and T are the 
infinites of space and time, respectively, in level a. For any time-space point 
{x, t) in level a, the monad of this point is expressed by {x, t) and the dimen- 
sions of the monad are (ei, £2, £3, Here ei,e2,Ss,et are infinitesimals of 
space and time, respectively, in level a. Moreover, assume that: 

= £2i^2 = l, £3^^3 = 1, StT^l (1) 

The space-distance between any two points 

Al = ^J{Ax,y + (AX2)2 + (AX3)3 

in level a, and 

Al' = ^J{Ax[f + (A4)2 + (A4)2 

in level p. The actual length (i.e., the length observed from the angle of level 
a) of Al' is: 
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II e-i — 82 — — e the actual length of A/' is: 



(2) 



(3) 



Similarly, the actual value (i.e., the value observed from level t) of At' is: 



And the actual lengths of the whole space and time axes in level (3 are: 



They are just the dimensions of a monad in level a. 

Now the limit of tending to zero (Axj — > 0) in level a is not tending to 
absolute zero but infinitesmal e. Here Axj can not go into the inner of the 
zero monad. In the inner of the zero monad in level a exists Ax^ instead of 
Axj. In other words, the process of tending to zero does not mean tending to 
absolute zero but the dimension of a monad in the same level. The standard 
part of any number u) is denoted by Stuj. To show clearly the meanings of 
Axj — > and At — > mentioned above, we do not use the notation Axj — > 
and At — > 0. Instead, we use 



(4) 



Lis\ — £i, ^2^2 — ^2-1 L^e\ — £3, Tel ~ 



(5) 



Sti^t 







(6) 



I.e., 




lim At = £t (7) 
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In other words, the definitions of {Axi £i) and (At £t) are, respectively, 
{StAxi —>■ 0) and {St At —>■ 0). Similarly, the limit of tending to infinite does 
not mean tending to absolute infinite but the nonstandard number L, i.e., 

lim Aa;i = Li, lim Ax2 — L2, lim Ax^ — L3, lim At — T 

St Axi-*oo St Ax2-*oo StAis— >oo 5tAt— >oo 

(8) 

Moreover, in standard function is expressed by f{xi,X2, X3, t) which 

represents the function value at the standard point {xi,X2,X3,t). However, in 
the nonstandard analysis, a point {xi, X2, x^, t) is a monad within which there 
are still infinite nonstandard points {x[,X2,x'^,t'). Then a function should 
be expressed by f{xi,X2,X3,t,x[,X2,x'^,t') which is the function value at the 
nonstandard point ( t') contained in the monad {xi,X2-,x^,t). Still 

we define the partial derivatives in form: 



dx'i stAx',^0 Ax', dt' stAt'^o At' ^ ' 



a/^ ^.^ f{x, + Ax,) - fix,) ^.^ fit + At) - fit) 

dxi stAxi-^o Axi ' dt stAt-^o At 



(i = l,2,3) 

Definition (9)-(10) is only a definition without the ordinary meaning of deriva- 
tive in the standard case. The difference of fix, + Axi) — fix,) in (10) is the 

^jbbrGVlQjtlOIl of f (^-^i ~i~ ^^-^i^ 5 '^kj ^i "^i^ j ' "^h^ ^ ) f^'^ii "^ji "^h^ ^3 "^ii ^ )' 

the increment of the function value between two monads of the points {xi+Axi) 
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and (xj). The distance between the points 

lim + Axi,^.) 

StAxi—^O 

and 

is infinite from the angle of level j3. Therefore, ^ and ^ in (10) have the 
meanings of an average. If one insists on the concept of ordinary derivative in 

standard case, ^ and ^ in (10) have no sense. 

2 Nonstandard Model Of Turbulence 

On the basis of the above-mentioned concepts of nonstandard analysis a 
nonstandard description of turbulent field can be presented as follows. 

Assumption 1: Global turbulent field is composed of standard points, 
and every standard point is yet a monad. Each monad possesses 
the internal structure, namely a monad is also composed of infinite 
nonstandard points (so called interior points of the monad). 

It is well known that a space is formed from the points. In fluid mechanics, 
fluid mechanists always take the fluid element (fluid particle) as point. In fact, 
a fluid particle has the volume not to be absolute zero but to be a micro volume, 
which is taken as infinitesimal from macroscopic viewpoint and infinite from 
microscope. People generally think of the whole fluid particle as uniform. 
Namely the fluid particle is taken to be equal to a point in reality. Now, in 
the case of turbulence, the motion of fluid varies so fast that a fluid particle 
can not be uniform as a whole but becomes a monad. Prom now on in this 
paper, the meanings of a fluid particle, a monad on a whole and a standard 

6 



point are identical with each other and so are the meanings of a fluid particle 
in the lower level, an interior point of a monad and a nonstandard point. 

Since the point of the turbulent field becomes a monad with internal struc- 
ture, the motion features of different interior points of the monad are different. 
So it is obvious that there is a flow or flow fleld in the inner of the monad. 
Hence, a fluid particle (monad) is made up of nonstandard points. The vol- 
ume of these nonstandard points is not absolute zero, and these nonstandard 
points, in fact, are the fluid particles in the lower level. Now the particles 
in the lower level are thought to be uniform. They are points (nonstandard 
points) actually. Moreover, a lot of fluid molecules are included still in the 
fluid particle of the lower level. 

Therefore, two kinds of flow fields exist in turbulence. One is the global 
turbulent field and the other is the fiow field in a monad ( called monad field 
in this paper). They are fiow fields in two different levels. Now the global 
turbulent field is not composed of the points in which there are no structures, 
but of monad fields. So there are three levels in turbulent field: molecular 
movement, monad field and global field. 

Assumption 2: The flows in monad fields are controlled by the Navier- 
Stokes equations.. 

Then the motion-equations of the fiows in monad fields can be written as 
follows (for incompressible fiuid, unsteadiness and three dimensions): 



dUi muj _ 1 doij 



Here f/j is velocity component in i-direction, p density of fluid, Ujj stress 
tensor and aij = aji. The independent variables of these functions in (11)- 
(12) are {xi, X2, X3,t, x[, X2, x'^,t'). et, 81,62, 83 are inflnitesimals (i.e., the lin- 
ear dimensions of a monad) of time t, space Xi,X2 and X3, respectively. Let 
et — El — 82 — S3 — e here and later, then the governing equations becomes: 

dUj dUjUj _ Idajj 

dt' dx'j ~ p dx'j ^ ' 

The governing equations (11)- (12) show the internal structure of a monad. 

For actual turbulence the possible linear dimensions of flow flelds in two 
levels could be estimated roughly as follows: 

Global turbulent held: inflnite L — ~ 10^-^cm 

inflnitesimal e ~ 10~^'^ — lQ~^cm 

(linear dimension of a monad) 
Monad fleld: inflnite e ~ 10"^'^ — 10~^cm 

inflnitesimal ~ 10~^'^ — W~^cm 

(linear dimension of the particle of lower level) 
According to this estimation, one space monad of the turbulent fleld con- 
tains approximately 10^ — 10^*^ nonstandard points (the particles in the lower 
level) . If the fluid is gas, there are approximately 10^ — 10^ gas molecules in one 
nonstandard point of a monad fleld. This number of gas molecules is enough 
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for stable mean values in the statistical average. 
Assumption 3: Turbulent field is continuous. 

Here the meanings of continuity are: firstly the global turbulent field is 
composed continuously of standard points (i.e., monad fields); secondly the 
flows in the monad fields are continuous; and thirdly the physical quantities 
on the interface between two infinitely close monads are continuous as well. 
The monads of 

lim {xi + Axi , Xj ,Xk,t) 

StAxi—^O 

and 

(Xj, Xj, X/j, t) 

are taken as two infinitely close monads in space. And so are in time. The 
mathematic expression of the third is 

lim U(x + Ax,0) = U(x,L), lim U(t + At,0) = U(t,T) (15) 

St Ax-^O ^ ' ^ ' StAt-*Q ^ ' ^ ' ^ ' 

There is an even more important meaning of the continuity. That will be 
stated in Assumption 6. 

3 Uncertainty Of Measurement 

To make a physical measurement, some time-space point (xi, X2,xs, t) has 
to be assigned first, at which the measurement is carried out. The results of 
the measurement are taken as the measuring data of this point {xi,X2,Xs,t). 
According to the standpoint of nonstandard analysis, the point {xi,X2,X3,t) 
is not an absolute point but a monad, in which there are infinite nonstandard 
points. Generally speaking, the motion features of different nonstandard points 
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are different. So tfie meaning of saying "measuring at point (xi, 0:2, a^s, t)" 
becomes ambiguous. Obviously, if tliere are fields in two different levels in the 
system being studied, there must exist the new meaning of a measurement at 
any point of the field in higher level. 

Assumption 4-' When a measurement at any point (monad) {xi,X2, xs, t) 
in a physical field is taken, the operation of the measurement will 
act randomly on one interior point (nonstandard point) of the point 

{xi,X2,X3,t). 

This assumption states clearly that one measurement operation at some 
standard point (monad) in a field will always take place at an interior point 
of such a monad. But the measurer can not determine which interior point is 
acted on by the measurement operation. What actually happens is random. 
The measurer knows only the standard point (monad) in which the interior 
point acted on by one measurement operation is contained. It is surprising that 
when a measurement is carried out there is no way to determine the exact po- 
sition of the measured object. That is the uncertainty of measurement. At 
first sight it seems that the uncertainty of measurement stems from the rough 
measuring skill and that the circumstance conditions of the measurement are 
not exactly controllable etc.. After further thinking, it is comprehended that 
the uncertainty rises essentially from limited ability and means of the descrip- 
tion of the physical world in hierarchical structures. Apparently, in their quest 
of the laws of the universe, people need not (and can not in fact) describe all 
motion states of every molecule and electron in the universe. Conversely, even 
though an exact and perfect description of the motion states of every molecule 
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and electron in the universe is given, the laws of the motion of the whole uni- 
verse are still not understood. It is clearly that because the physical world has 
a hierarchical structure, therefore any description and measurement of phys- 
ical phenomena should be related with the hierarchical structure. When the 
effect of hierarchical structure can not be ignored, the uncertainty of measure- 
ment occurs. The measurement and observation of turbulence can be taken 
for examples of this case. 

An operation of measurement will act randomly on various interior points 
in a monad. By what probability does the measurement operation act on each 
interior point of a monad? In turbulence it is thought that one operation 
of measurement will take place in equiprobability at each interior point of a 
monad. 

Assumption 5: When a measurement at any point (monad) of a tur- 
bulent field is made, the operation of the measurement will act in 
equiprobability on various interior points of the monad. This as- 
sumption is called the equiprobability assumption. 

So there are two kinds of averages. One is that one measurement is just 
taking average over a large number of molecules contained in some interior 
point of a monad, while the other average over the motions of all interior points 
in a monad has to be taken in order to get the mean values of the physical 
quantities over the standard point (monad) of the turbulent field. According to 
the equiprobability assumption, the second average formula can be given. Let 
U {x, t, x\ t') represent some physical quantity at an interior point (x', t') of the 
monad {x,t) and U{x,t) express the average of U{x,t,x',t') over all interior 
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points of the monad. Then it follows that 



U(x,t) = -f dt'—r[ dx[ [ dx'2 I dx'M(x,t,x',t') (16) 
T Jo L'^ Jo Jo Jo 

Where the monad {x,t) = {xi,X2,X3,t), the interior point of the monad 
{x',t') = {x[,X2,x'^,t'), and T, L are infinites of time and space respectively. 
Let T = Li = L2 = L:^ = L here and later. Obviously, U is only a function of 
the standard point {x,t). 

If the monad fields of a global flow field are uniform, there is no difference 
between saying "one measurement acting on the interior point of a monad" 
and "one measurement acting on the point (monad)". Then the global field is 
laminar. Exactly,the essential feature of turbulence is that its particle (stan- 
dard point) has internal structure, i.e., internal flow field (monad field). That 
could be taken as the definition of turbulence. 

Because of the uncertainty of measurement, one measurement at the point 
{x, t) will act randomly on any interior point. Therefore, the results of many 
measurements carried out at the same point (monad) under constant circum- 
stance conditions will indicate randomly the motion states of various interior 
points of the monad. Surely fluctuation of data of the measurements must 
appear, provided that there are considerable different states among various 
interior points of the monad. This is the mechanism of fluctuation appearing 
in a turbulent field. 

Thus, a turbulent field is composed of monads with interior structure 
(monad field). The flows in the monad fields are controlled by the govern- 
ing equations of (13)-(14). Still continuity of turbulent field exists. Therefore, 
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turbulence is also regular flow. The appearance of the fluctuation of a turbu- 
lent flow stems from the uncertainty of measurement and observation in the 
turbulent fleld. 

To keep the concepts mentioned above in mind, the order of magnitude 
of fluctuation in a turbulent field can be estimated roughly. Let Ax' be finite 
length in the coordinate system of a monad, U ~ 0(1) any physical quantity 
and its fluctuation u ~ U{x' + Ax') - U{x') - 0Aa;'. The actual rate of 
change (i.e., observed from the angle of the global fleld level) of function U 
is In one case |^ ~ 0(6^), u ~ 0(e:^), no obvious fluctuation occurs. 

In another case ^ ~ 0(1), and ~ O(^) showing that there is excess 
rate of shear strain in the fleld. This will lead to instabihty in the flow fleld. 
So the case should be excluded. Then in the third case ^ ~ 0(e) only, and 
~ O(^) with fluctuation u ~ 0(e). So there is notable fluctuation in 
order of magnitude u ~ 0(e) in turbulence provided that there exists any set 
of interior points in a monad, and on the set |^ ~ 0(e). 

4 Conservation Equations 

Generally speaking, in fluid mechanics a small volume Axi Aa;2 Axa is taken 
and the flux of conservative physical quantity through the boundary planes of 
the volume is computed. Then let Axi 0, Ax2 0, Axs — > 0, the small 
volume will tend to the point (xi, X2, X3). Thus, conservation equations at 
the point {xi,X2,Xs) are obtained. The boundary planes of the small volume 
are absolute planes without thickness. However, in nonstandard case a small 
volume is still taken. By Use of the same method as in standard case but 
let StAxi 0, St Ax2 —>■ 0, StAx^ — > 0, the small volume Axi AX2AX3 will 
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become the monad of {xi, X2, x^). Now the boundary planes of the small 
volume are not absolute planes but very thin layers with thickness of e (the 
linear dimension of the monad) . Also conservation equations at this monad are 
obtained. Surely the equations are those of mean quantities over the monad 
in which infinite interior points are contained. 

Now these conservation equations are given in the following. 

Let " ~ " express the average operation over {x[, X2,x'^, t'), " ~* " over x[, 
" — " over t', etc., in a monad. 

Through the pair of the planes perpendicular to Xi axis the mass inflow 
per unit time is: 



- AX2A.X3 



(xi)Aa;2Aa;3 



the mass outflow 



{xi + Aa:i)Aa;2Aa;3 



the net inflow 



dxi 



d —I 



Ax2Aa;3 



Aa;iAa;2Aa;3 



Through the other two pairs of the planes perpendicular to the two other axes 
the net mass inflow per unit time are, respectively. 



U2 AX1AX2AX3 



dx 



and 
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-C/3 AX1AX2AX3 



Because of the law of mass conservation, it follows that 

„ '/^X2^X3 ^ ;Aa;iAx3 „ -Aa;iAa;2 

axi 0x2 ox^ 

Here the average of " ~ Aa;iAa;2" is taken over the plane of Aa;iAa;2, which 

tends to " ~ 12" when St/S.xi — > 0, StAx2 — > 0. And so are the others. Let 

StAxi 0, 6'tAa;2 — > 0, StAx^ — > 0, then it is obtained that 



d ^123 d ^123 (9 —123 



5xi dx2 dx 

i.e., 

atTi Sf/a dU^ ^ 

Similarly, the equations of momentum conservation can be obtained too. 
Through the pair of the planes perpendicular to Xi axis, the net inflow of 
the momentum in xi direction per unit time is 



pUiUi {xi)Ax2Ax3— pUiUi (xi + Axi) Ax2Aa;3 

— ■ Aa;2Aa;3 



d 



= -p^ — UiUi Aa;iAa;2Aa:3 

0x1 

Through the other two pairs of the planes perpendicular to the other two axes, 
the net inflows of the momentum in xi direction per unit time are 



-p-?^ — U2V\ AxxAx2Ax-i 
0x2 
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and 



0x3 

respectively. 

The force in xi direction exerting on the pair of the planes perpendicular 

to Xi axis is: 



-ou {xi)Ax2/^xs + CTir (xi + Axi)Ax2/^X3 

Q ^;^^A.X2Ax3 

-Oil Aa;iAa;2Aa;3 



dxi 

The forces in xi direction exerting on the other two pairs of the planes per- 
pendicular to the other two axes are 



^Xi^X2^X2, 



and 



Q j^gAa;iAa;2 

-oil AX1AX2AX3 



8X2, 

respectively. 

Then the increments of the momentum in xi direction in the small volume 
per unit time are 



~_AxiAx2Axz ^Aa:iAa;2Ax3 AxiAxtAxs 

p{U, {t + At) - C/i {t)} 

Q _2:^AxiAa;2Aa;3 

= p-^Ui Aa;iAx2Ax3 
Because of the law of momentum conservation, it follows that 
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p—Ui +Pt^UiUi +Pt^U2Ui +P^U3Ui 

at oxi 0x2 0x3 



Q ;^^AxiAxi Q ^^^AxxAxo, Q ^^^Ax\Axi 

= ^ — ^ + ^ — ^ + ^ — ^ 

OX\ 0X2 0X3 

Let StAxi 0, StAx2 0, StAx3 0, so " ~ Axi AxaAxg" tend to " ~ 
123" , etc. Then it follows that 



a ^123 d ~ d ~ a ~ 

-t/l +7^t/lC/l +l^f/2f/l +^t^3C/l 
or 0X1 (7X2 0x3 



1 a ~l23 1 d —2^3 1 5 —312 

— c^ii + -t; — cr2i + -t; — ctsi 
/9 0x1 p 0x2 p 0x3 



I.e., 



O r\ o o 

ar axi ax2 

= -(7^ — 0-11 + 7- — a2i + ^ — 0-31) 

p 0x1 0x2 0x3 

Similarly, the momentum equations in X2 and X3 directions can be written 
too. The general form of momentum equations is: 

dUi dUiUj Idd^j 
dt dxj p dxj 

The equations (17) — (18) are about the mean quantities over a monad and 
called the conservation equations. Here no special mean-methods are assigned. 
But the mean " ~ " will always be taken as the mean- method showed in (16) 
in this paper. 
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5 Fundamental Equations of Turbulent Flows 

Still instantaneous quantities Ui and P are decomposed into two parts such 
that 

Ui^lJi + Ui, P^P + p, Ui^O, p=0 (19) 

Here the quantities with " ~ " are called mean quantities; in this regard, 
we will adopt the mean- method of (16) i.e., 

[/ = - / dt'— / dx[ / dx', / dx'.U, U^- dt'U, U'^- dxp 
1 Jo Jo Jo Jo 1 Jo L Jo 

(20) 

= ir^ £ "< £ = ^ £ "< £ £ ^^'^ 

The relations of (19) are the definitions of Ui and -p. The ui and p have the 
same the order of magnitude as the turbulent fluctuation mentioned in Section 
3. We will call the Ui and p in (19) as fluctuant quantities. To note that the 
f/j and P are mean quantities over a monad, and Ui and p are meaningful to 
a monad (or in a monad). They have different meanings with those in the 
existent theory of turbulence. 

Before derivation of the fundamental equations of turbulent flows, it is 
necessary to present the assumption of close property between two monads 
when they are infinitely close to each other, and analyze the mutual relation 
of the two monads in motion characteristics. 

Assumption 6: In both the value and structure of function, physical 
function, defined on the interior points of the monads of a turbulent 
field, is very close between two monads, when these two monads are 
infinitely close to each other. 
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Considering any physical quantity Ui ~ 0(1), its fluctuation Ui, the math- 
ematical expressions of this assumption are as follows: 

(a) dui - 0(£2), (6) d{uiUj - ti^j) ~ 0(£^) (21) 
(0 < < L, < t' < T) 
Here the operator 6 is defined as 

S^U= lim [U{x,+Ax„x'^-U{x,,x% 6tU = lim [U{t+ At, t')-Uit,t' )] 

olAa;,— >U ocAt— >U 

(22) 

It can be got from (a) in (21) that 

S \ - — 
e^dx'i dxi 

Where is the actual rate of change or figure slope of the function U at 
the nonstandard point x[, and |^ is the mean value of over (see(25)). 
Clearly the condition (23) and (a) in (21) show that the figure shape of the 
function U has very little difference between two infinitely close monads. Sim- 
ilarly the condition (6) in (21) shows that the figure shape of the function UiUj 
has very little difference between two infinitely close monads. 

Assumption 6 could be called "Continuous Assumption" . It is well known 
that with a few exceptions, physical quantities vary usually continuously with 
time and space. In the case of standard analysis, a function is defined on 
the standard points. It is characterized only by its value. The continuity of 
a physical quantity only requires that the function representing the physical 
quantity be continuous in mathematics. On the other hand, in the nonstandard 
case a point becomes a monad and a function is defined at the interior point 
of the monad. There exists a function-structure on the whole monad. So 
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0(£^) 



(23) 



the function is characterized by not only its value but also its structure on 
the monad. The continuity of a physical function will mean that in both the 
values and figure-shape of the function , the variation of the physical function 
with monads is continuous. Thus, there is very little difference between two 
infinitely close monads in both value and structure of the physical function. In 
other words, the Assumption 6 is the nonstandard expansion of the physical 
function continuity in the case of standard analysis. Moreover, the Assumption 
6 is reasonable and natural, because the circumstance conditions imposed on 
two infinitely close monads are very close to each other. 

In a monad, we consider that stress tensor for Newtonian fluid is 



Using (20) and(17),(15) (Assumption 3) we can give the following computation: 




(24) 




(25) 



Then 






By use of (a) in (21) it follows that 




(26) 
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To omit them, it is obtained that 

da^ _ _dP_ d^Uj 
dxj dxi ^ dxj 

Then the equations (18) can be written as 

dUi dUiU. IdP „2fr ^2 d d 

V + = + V2 = — — 28 

Ot OXj p Xi OXj OXj 

(17) is written as 

Here ^ — ^ the kinematic viscosity. 

Note that now the mean method in the equations (28) and (29) is defined 
by (20). Then by use of the relation (19), the decomposition of (29) and (28) 
can be done as 



KJJU't KJJU't 



dt dxj dt dxj 

= --1^ + + ^^"U, - (31) 
p OXi P dXi 

Here Ui,p and (UiUj) fi are the fiuctuations of the velocity, pressure and (UiUj) 
respectively. The independent variables of the functions in the equations (30)- 
(31) are {xi, X2, X3,t, x[, X2, x'^jt'). Based on Assumption 6 the order of mag- 
nitude of the terms concerned with fluctuation quantities is ~ 0(e), but the 
order of magnitude of the terms concerned with instantaneous quantities is 
~ 0(1). So when (30)-(31) are split into two parts having different order of 
magnitude, the following can be obtained: 
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Equations, in which terms are in the order of magnitude ~ 0(1), of instan- 
taneous quantities are: 

■air' '^'^ 

f + fa^-i|^ + .VV. (33) 

at oxj p oxi 

Equations, in which terms are in the order of magnitude ~ 0(£), of fluctu- 
ation quantities are: 

air' 



dui d{UiUj)fi 1 dp „o 
at oxj p oxi 



Where {UiUj)fi in (35) is: 



{UiUj)fi = UiUj-UiUj 

— UiUj + UjUi + UiUj — UjUj 

Substitution of the relations into (35) produces 

duj ^ ^dUj ^ ~ duj ^ dujUj du^j ^ 1 ^ ^.^2^ /ggN 
dt ^ dxj ^ dxj dxj dxj pdxi * 

Equations (28) (29), (32) (33) and (34) (35) are the governing equations about 
mean, instantaneous and fluctuant quantities respectively. They can be called 
as the fundamental equations of turbulent flows. 

It is found easily that these equations have the same form as those in the 
existent theory of turbulent flows. But there are essentially differences between 
the two. These differences are: 
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1. Different metfiods to get tfiem. In tfie existent tfieory of turbulence the 
mean equations are obtained by tlie average of tlie instantaneous equations 
(i.e., the Navier-Stokes equations). Then the mean equations are subtracted 
from the instantaneous equations to get the fluctuation equations. But in this 
paper, the mean equations will FIRST be obtained in view of the physical 
conservation laws, and then by virtue of the close property assumption, the 
instantaneous and fluctuant equations are obtained respectively. 

2. The terms in the equations of the existent theory are the limits, such as 

^.^ fix + Ax) - fix) fit + At) - fit) 

Ax^O Ax At^O At 

But the terms in the equations of the nonstandard theory are the partial 
derivatives defined in (10). The equations in two theories are conceptually 
different. The former is based on, so called, the frame oi 6 — e, the latter is 
out of the frame of S — e. Therefore, the equation in nonstandard picture of 
turbulence is the new kind of equation. 

3. The physical quantities in the equations of the existent theory, such 
as velocity, pressure etc., are all functions of the standard point ix,t). How- 
ever, those in the nonstandard case are all functions of the point ix,t,x',t'), 
actually functions of the nonstandard point ix',t') with ix,t) showing which 
monad contains it. It is reasonable that the physical quantities are taken as 
the functions of the standard point (x, t) in laminar flow, but not in turbulent 
flows, because a standard point becomes a monad in which there is an interior 
structure (i.e., the monad flow field) in the case of turbulence. In fact, there is 
no definition of instantaneous and fluctuant quantities of a whole monad, but 
there is the definition of instantaneous and fluctuant quantities of the interior 
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points of the monad. So for the global field,the Navier-Stokes equations in 
existent meanings do not exist in principle. Therefore, though (32)-(33) are 
the same as the Navier-Stokes equations in forms, they are different. The fun- 
damental equations obtained in this paper arc governing equations about the 
physical quantities defined at those interior points of various monads. These 
interior points have the same nonstandard coordinates {x',t'). But the fun- 
damental equations in the existent turbulent theory are governing equations 
about the physical quantities defined at the standard points of the global flow 
field. There is an essential difference between them. 

4. The average in the existent theory usually means the ensemble aver- 
age, but the distribution function can not be given. So the closure problem 
continues to exist. Otherwise, it is definite in the nonstandard model to take 
average over all interior points of a monad by the equiprobability assumption. 

5. It should be pointed out that the function values at some standard 
points are given for the initial and boundary conditions of the fundamental 
equations in the existent theory. However, only the function values at some 
nonstandard points can be taken as the initial and boundary conditions for the 
instantaneous and fiuctuation quantities in the fundamental equations in the 
nonstandard case. Because a turbulent field always starts and develops from 
walls or laminar fiows and evolves from laminar fiows or static fiuids, the func- 
tion values on walls or laminar fiows can be taken as the boundary conditions, 
and the conditions of static fiuids or laminar fiows can be taken as the initial 
conditions. No slip condition is still kept. Thus, the initial and boundary con- 
ditions of the fundamental equations obtained from the nonstandard analysis 
can be given in principle. 
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6 On Closure Problem 



Substitution of the relation Ui = Ui + Ui into (28) produces 



dUi dUiUj duiUj 




(37) 



dt dxj dxj 



This is the same in form as famous Reynolds equation. It is not closure, and 
neither is equation (36) in which there is term UiUj. The closure problem 
stemming from nonlinearity is very difficult in turbulence research. Now in 
the nonstandard theory, the turbulence is composed of monad fields, and the 
point-average is adopted. Then by the close property between two infinitely 
close monads, we can give the reasonable closure methods as follows. 

By the relation (21) (Assumption 6) and the order of magnitude of fluctu- 
ation Ui ~ 0(e), it follows that 



5{uiUj) — Ui5uj + UjSui ~ 0(£^) 




duiUj 



O(e^) 



and 



5{uiUj — UiUj) ~ 



also, it follows that 



d{uiUj — UiUj) 



(38) 



dxj 



There are three choices for the closure of equations: 
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(A) Choice one: 

The term ~ 0(£^) is neglected from equations (37). The closed equa- 
tions of the mean quantities will be obtained as follows. 

(39) 



dxi 

dUj dUjUj _ _ldP_ 
dt dxj pdxi ' 



(B) Choice two: 
To note that 



d{UiUj — UiUj) 3 



dxj 



Q{s') (38) 



and the other terms in equation (36) have the order of magnitude of ~ 0(e). 

Hence, the term (38) can be neglected from (36). Then the equation (36) 
becomes 

^ +u— + IJ — = --— + uV^u- (41) 
dt ^ dxj ^ dxj pdxi * 

With Ui — Ui — Ui, this equation changes to 

dui dui dUi dui 1 dp 2 /.o^ 

-wr + Uj^ + uj- 2uj—^ — — + iyVui (42) 

at oxj oxj oxj poxi 

Under certain initial and boundary conditions, the equations (32)-(33) are 

solved and Ui are obtained. Then the equations (42) with (34) becomes also 

closed, the fluctuations Ui and p can be obtained after the equations (34) (42) 

are solved. Finally the mean quantities Ui = Ui — Ui and P = P — p. 



(C) Choice three: 
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Hence, the term (38) can be neglected from (36) and (37). Then it follows 
that 

■air' (^^) 



+ ^^ + ^^ = —T— + i^^^Ui (44) 



dt dxj dxj p dx 

dui 



dui ~dui dUi 1 dp 2 

'dt^^^d^-^ ""^d^- ^~~pd^-^^^ ^^^^ 

Obviously, the equations (43)-(46) are closed. We can obtain Ui,P,Ui,p, 
from solving these equations. 

Surely ii Ui — 0,p — 0, i.e., no fluctuations occur, the equations (28) (29) 
and (32) (33) are reduced to ordinary Navier-Stokes equations and the field 
becomes a laminar field. 

7 Some Remarks on the Numerical Calculations 

To take the ^ for example, in the standard case we have 

( ^\ - ^(^1^^2,3:3,^ + At) - U{xi,X2,X3,t) 

but in the nonstandard case by the definition of ^ in this paper, it follows 

that 

(—\ - lim ^i^^^^2,X3,t + At,x[,x'2,x'^,t') - U{xi,X2,Xs,t,x[,x'2,x'^,t') 
\dt J ~ stM^o At 

(47) 

In form, there is no difference between the two cases when the discretization 
of time-space is taken for the numerical computation. However within a grid, 
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obtained from the discretization of j ^, there is not structure (i.e., the whole 
grid should be uniform). But the interior structure is permitted through the 
grid obtained from the discretization of (^^^ ^- This diference between two 
cases is very important. 

Surely in the nonstandard case, U {xi, X2, Xs, t+At, x[,X2, Xg, t') and U {xi,X2, X3, t, x'l, x'2, x'^, t') 
are the values of physical quantity U at the interior point ( t') of 

two different monads in time, that is monad of {t + At) and monad of {t). 
Though the fundamental equations work for any interior point {x[,X2,x'^,t') 
of every monad, it is better that let x[ = X2 = x'^ = t' = because the 
point {xi, X2, X3, t, 0, 0, 0, 0) is located on the boundary of the time-space of a 
turbulent field, when {xi,X2,X3,t) is the boundary monad of the field. 

Moreover, U {xi,X2, X3,t + At, 0, 0, 0, 0) and U {xi,X2, X3,t, 0, 0, 0, 0) are hy- 
perreal numbers. In the numerical calculations, their standard parts will be 
taken. 

Therefore, after the discretization, ^ becomes: 

StU{xi,X2, X3, t + At, 0, 0, 0, 0) - StU{xi, X2, xg, t, 0, 0, 0, 0) 

Xt 

For simplicity it can also be written as 

U{xi,X2, xs, t + At, 0, 0, 0,0) -U {xi,X2, X3, t, 0, 0, 0, 0) 

At ^ ^ 

provided it is kept in mind that the U in (48) represents the standard parts of 
them. 

There are similar discussions for the space partial derivertives of physical 
quantities. 

After the U{xi,X2,X3,t,0,0,0,0) and u{xi,X2, X3,t,0,0,0,0), the instan- 
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taneous and fluctuant quantities respectively, are obtained, then 

U{xi,X2,X3,t) = U{xi,X2,X3,t,0,0,0,0) - u{xi,X2,X3,t,0,0,0,0) (49) 

That is the mean value over the small volume of AxiAx2A.X3At, but it is taken 
as the mean value over the point {xi,X2, X3,t). 

Moreover, the mean value U obtained by the methods mentioned above, 
is the mean over the point (monad) and still random oscillatory. The random 
oscillation of the mean values U ought to stem from the unavoidable random 
disturbances in real turbulent flelds. If the average of these mean values over 
a finite time period(or a finite space range) is taken once again, the final 
average results could be compared with the measured average values of physical 
quantities over corresponding time period(or space range). 

8 Conclusions 

In a view of the concepts of the nonstandard analysis and the fact that 
the physical world has a hierarchical structure, one nonstandard picture on 
turbulent fiows has been presented in this paper. The key points of this picture 

are as follows: 

1. There exist two kinds of fields in different levels. One is the global 
turbulent field composed of the standard points (monads), and the other is 
the monad field. And as a whole a monad is not already uniform. Now the 
interior structure occurs in a monad; namely, different interior points of the 
monad possess different characteristics of motion. 

2. The fiows in a monad field are governed by the Navier-Stokes equations. 

3. If there are two kinds of fields in different levels, one operation of 
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measurement at any point (actually a monad) of the field in a higher level 
will act randomly on one interior point of the monad field. The measurer can 
determine the monad, rather than its interior point on which the measurement 
operation acts. 

The position of measured object can not be determined exactly, but over a 
range. Evidently there exists some uncertainty of measurement. That uncer- 
tainty stems from the description and observation of the physical phenomena 
in a hierarchical structure. 

4. The flows in a turbulent field are also continuous and regular, while 
they are seen intuitively in disorderly and irregular form. Fundamentally, 
the disorder, irregularity and fiuctuation occur because of the uncertainty of 
measurement. Hence, we can probably say that the disorder of turbulence 
stems from two sources: the first is the fiuctuation because of the uncertainty 
of measurement. This kind of fiuctuation is thought to be "real turbulent 
fiuctuation". The second is the violent unsteadiness and unhomogeneity of 
the field. The values of Ui and P vary fast with time and space. It is due 
to the flow-instability as the Reynolds number is very high. In this case, the 
small random disturbances will be amplifled. 

5. There are two kinds of averages. The values of the measurement acting 
on one interior point of a monad, merely stand for the mean results over a lot 
of fluid molecules contained in this interior point. The other is the average 
over all interior points in a monad. The second average results indicate the 
mean values of the physical quantities over the monad (the standard point) in 
the global fleld. By the equiprobability assumption, the mean formula of the 
second can be given in (20). 
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6. As the expansion of the fact that the physical functions are continuous 
in the standard case, the assumption of close property between two infinitely 
close monads shows that the variation of interior structure as well as values 
of the physical function U{x,t,x' ,t') with {x,t) (monad) is continuous too. 
In other words, there exists only very little difference in characteristics of 
the function U{x,t,x',t') between two infinitely close monads. According to 
this assumption and the conservation laws, the governing equations of mean, 
instantaneous and fluctuant quantities in turbulence (i.e., the fundamental 
equations) are obtained, and the reasonable closure methods are given too. 
It should be noted that the point-average is used, instead of the Reynolds 
average, therefore the closure problem is easily overcome in the nonstandard 
case. 

7. The new kind of equations of turbulent motion is presented in this paper. 
These equations are out of the frame oi S — e and valid for the points, which 
are permitted to be not uniform. These points, in fact, are monads. 
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